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METHOD OF IDENTIFYING STRATEGIES FOR TREATMENT OR 
PREVENTION OF VENTRICULAR FIBRILLATION AND VENTRICULAR 

TACHYCARDIA 

[0001] This application is a national stage application vmder 35 U.S.C. § 371 

from PCX Application No. PCT/US2005/002848, filed February 2, 2005, which 
claims priority benefit of U.S. Provisional Patent Application Serial No. 60/541,531, 
filed February 3, 2004, which is hereby incorjDorated by reference in its entirety. 
[0002] This appHcation was made, at least in part, with funding received from 

the National Institutes of Health under Contract No. HL62543. The U.S. government 
may r e tain has certain rights in this invention. 

FIELD OF THE INVENTION 

[0003] The present invention relates to a method of identifying strategies for 

the treatment or prevention of ventricular fibrillation and ventricular tachycardia. 

BACKGROUND OF THE INVENTION 

[0004] Catastrophic heart rhythm disorders are among the leading causes of 

death in the United States, The most dangerous of these arrhythmias is ventricular 
fibrillation, a disturbance in which disordered wave propagation causes a fatal 
disruption of the synchronous contraction of the ventricle. Although the exact 
mechanism for fibrillation is still being debated, one theory proposes that fibrillation 
is a state of spatiotemporal chaos consisting of the perpetual nucleation and 
disintegration of spiral waves (Chen et al. , "Mechanism of Ventricular Vulnerability 
to Single Premature Stimuli in Open-Chest Dogs," Circ Res 62:1991-209 (1988) and 
Witkowski et al,, ''Spatiotemporal Evolution of Ventricular Fibrillation," Nature 
392:78-82 (1998)), in association with a period doubling bifurcation of local electrical 
properties (Garfinkel et al., "Preventing Ventricular Fibrillation by Flattening Cardiac 
Restitution," Proc Natl Acad Sci USA 97:6061-6 (2000); Gilmour et al., "Electrical 
Restitution, Critical Mass, and the Riddle of Fibrillation," J Cardiovasc 
Electrophysiol 10:1087-9 (1999); Karma, A., "Electrical Altemans and Spiral Wave 
Breakup in Cardiac Tissue," Chaos 4:461-72 (1994); and Panfilov, A., "Spiral 
Breakup as a Model of Ventricular Fibration," Chaos 8:57-64 (1998)). Nucleation of 
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the initiating spiral wave pair is caused by local conduction block (wave break) 
secondary to spatial heterogeneity of refractoriness in the ventricle (Arce et al,, 
"Altemans and Higher-Order Rhythms in an Ionic Model of a Sheet of Ischemic 
Ventricular Muscle," Chaos 10:41 1-26 (2000); Chen et al, "Mechanism of 
Ventricular Vulnerability to Single Premature Stimuli in Open-Chest Dogs," Circ Res 
62:1991-209 (1988); Sampson et aL, "Simulation and Prediction of Functional Block 
in the Presence of Structural and Ionic Heterogeneity," Am J Physiol 28 1 :H2597-603 
(2001); Winfree, When Time Breaks Down (Princeton, NJ: Princeton University 
Press) (1987); and Witkowski et al,, "Spatiotemporal Evolution of Ventricular 
Fibrillation," Nature 392:78-82 (1998)). Until recently, spatial heterogeneity was 
thought to result solely from regional variations of intrinsic cellular electrical 
properties (Sampson et al., "Simulation and Prediction of Functional Block in the 
Presence of Structural and Ionic Heterogeneity," Am J Physiol 281 :H2597-603 

(2001) ; and Yan et al., "Characteristics and Distribution of M Cells in Arterially 
Perfused Canine Left Ventricular Wedge Preparations," Circulation 98:1921-7 
(1998)) or from stimulation at more than one spatial location (Qu et aL, "Mechanisms 
of Discordant Altemans and Induction of Reentry in Simulated Cardiac Tissue," 
Circulation 102:1664-70 (2000); Watanabe et aL, "Mechanisms for Discordant 
Altemans," y Car<^/ova5c Electrophysiol 12:196-206 (2001); and Winfree, When Time 
Breaks Down (Princeton, NJ: Princeton University Press) (1987)). However, it is now 
appreciated that purely dynamical heterogeneity can be sufficient to cause conduction 
block during single-site stimulation in both homogeneous one-dimensional models of 
canine heart tissue and in rapidly paced canine Purkinje fibres (Echebarria et al., 
"Instability and Spatiotemporal Dynamics of Altemans in Paced Cardiac Tissue," 
Phys Rev Lett 88:208101 (2002) and Fox et al, '^Spatiotemporal Transition to 
Conduction Block in Canine Ventricle," Circ Res, 90:289-96 (2002)). A similar 
mechanism has been shown to precipitate conduction block and spiral break-up in 
models of homogeneous two-dimensional tissue (Fenton et aL, "Multiple Mechanisms 
of Spiral Wave Breakup in a Model of Cardiac Electrical Activity," Chaos 12:852-92 

(2002) ). 

|0005| The period doubling bifurcation implicated in the transition to 

conduction block is manifest as alternans, a beat-to-beat long-short alternation in the 
duration of the cardiac action potential (Garfmkel et al., "Preventing Ventricular 



3 



Fibrillation by Flattening Cardiac Restitution," Proc Natl Acad Sci USA 97:6061-6 
(2000); Gilmour et al., "Electrical Restitution, Critical Mass, and the Riddle of 
¥\hm^i\on;' J Cardiovasc Electrophysiol 10:1087-9 (1999); Karma, A„ "Electrical 
Altemans and Spiral Wave Breakup in Cardiac Tissue," Chaos 4:461-72 (1994); 
Panfilov, A., "Spiral Breakup as a Model of Ventricular Fibration," Chaos 8:57-64 
(1998); and Watanabe et aL, "Mechanisms for Discordant Alternans," J Cardiovasc 
Electrophysiol 12:196-206 (2001)). Previous investigators have hypothesized that 
altemans can be accounted for by a simple uni-dimensional return map called the 
action potential duration restitution function (Chialvo et al., "Low Dimensional Chaos 
in Cardiac Tissue," Nature 343:653-7 (1990); Chialvo et al., "Non-Linear Dynamics 
of Cardiac Excitation and Impulse Propagation," Nature 330:749-52 (1987); Guevara 
et al,, "Electrical Altemans and Period Doubling Bifurcations," IEEE Comput Cardiol 
562:167-70 (1984); and Nolasco et al., "A Graphic Method for the Study of 
Alternation in Cardiac Action Potentials," J Appl Physiol 25:191-6 (1 968)). This 
hypothesis assumes the duration D of an action potential depends only on its 
preceding rest interval / through some function f (I) that is measured experimentally. 
If the D restitution function has a slope S:l , then a period doubling bifurcation occurs 
for some value of the stimulus period T, where T-D^L The velocity Fat which an 
action potential propagates can also be described by a restitution function, where V = 
c(I). 

[0006] It has also been shown previously that the combination of a steeply 

sloped action potential duration ("APD") restitution function and a monotonically 
increasing conduction velocity ('"CV") restitution function is sufficient to produce 
dynamical conduction block during sustained pacing at a short cycle length (Fox et al, 
"Conduction Block in One-dimensional heart Fibers," Phys. Rev. Lett. 89:198101 
(2002); Fox et al, ''Spatiotemporal Transition to Conduction Block in Canine 
Ventricle," Circ Res. 90:289-96 (2002)). This observation may provide a generic 
mechanism for wave break and the onset of ventricular tachycardia and fibrillation. 
However, it is unlikely that the conditions used to demonstrate this phenomenon 
experimentally apply to the clinical situation, where the induction of ventricular 
tachyarrhythmias typically is associated with the interruption of normal cardiac 
rhythm by only a few premature beats. A single premature beat is sufficient to cause 
spatial heterogeneity in the form of discordant altemans (Watanabe et al., 
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''Mechanisms for Discordant Ali^mans,'' J Cardiovasc Electrophysiol 12:196-206 
(2001)), but the conditions required for the development of conduction block in this 
setting have not been studied extensively. 

[0007] Other studies, building on earlier theoretical work by Krinsky, Winfree 

and colleagues (Krinsky et al, "Votices with Linear Cores in Mathematical Models of 
Excitable Media," PhysicaA 188:55-60 (1992) and Winfree, A., "Evolving 
Perspectives During 12 Years of Electrical Turbulence," Chaos 8:1-20 (1998)) and 
experiments by AUessie (Allessie et al., "Circus Movement in Rabbit Atrial Muscle as 
a Mechanism of Tachycardia. III. The 'Leading Circle' Concept: aNew^Modelof 
Circus Movement in Cardiac Tissue Without the Involvement of an Anatomical 
Obstacle," Circ Res 41:9-18 (1977)), have suggested that spiral wave re-entry could 
be the 'engine' that drives ventricular fibrillation ("VF") (Frazier et al., "Stimulus- 
Induced Critical Point. Mechanism for Electrical Initiation of Reentry in Normal 
Canine Myocardium," J C/f>7 Invest 83:1039-1052 (1989); Witkowski et aL, 
"Spatiotemporal Evolution of Ventricular Fibrillation," Nature 392:78-82 (1998); 
Weiss et al., "Chaos and the Transition to Ventricular Fibrillation: a New Approach 
to Antiarrhythmic Drug Evaluation," Circulation 99:2819-2826 (1999); Chen et al., 
"Mechanism of Ventricular Vulnerability to Single Premature Stimuli in Open-Chest 
Dogs," Circ Res 62:1 191-1209 (1988); Gilmour et al., "Electrical Restitution, Critical 
Mass, and the Riddle of Fibrillation," J Carc/Zova^c Electrophysiol 10:1087-1089 
(1999); and Pertsov et al., "Spiral Waves of Excitation Underlie Reentrant Activity in 
Isolated Cardiac Muscle," Circ Res 72:631-650 (1993)). Although there is substantial 
evidence that spiral wave re-entry contributes significantly to the induction and 
maintenance of VF, the exact mechanisms by which spiral waves sustain VF is 
currently being debated. 

[0008 1 The present invention is directed towards correctly identifying the 

mechanisms (or more likely, mechanisms) by which spiral waves cause VF and the 
development of pharmacological approaches to VF treatment and prevention. 



SUMMARY OF THE INVENTION 

[0009] The present invention relates to a method of evaluating the effect of a 

physiological condition on the occurrence of ventricular fibrillation. This method 
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involves providing a test system and initiating a ventricular fibrillation inducing 
sequence in the test system by interrupting normal sinus heart rhythm with premature 
electrical stimuli. Following the initiation of ventricular fibrillation a ventricular 
recovery sequence is initiated. The test system is subjected to a physiological 
condition before, during, and after initiating the ventricular fibrillation recovery 
sequence and the physiological conditions which affect ventricular fibrillation 
recovery are identified. 

[0010] Another aspect of the present invention relates to method of identifying 

treatment candidates as therapeutic strategies for the prevention of ventricular 
fibrillation. This method involves providing a test system and initiating a ventricular 

fibrillation inducing sequence in the test system by interrupting normal sinus heart 
rhythm with premature electrical stimuli before, during, and after administering the 
treatment candidate to the test system. Treatment candidates which prevent the 
initiation step fi*om inducing ventricular fibrillation are identified as therapeutic 
strategies for prevention of ventricular fibrillation. 

[001 1] A further aspect of the present invention relates to a method of 

identifying treatment candidates as therapeutic strategies for treating ventricular 
fibrillation. This method includes providing a test system and initiating a ventricular 
fibrillation inducing sequence in the test system by interrupting normal sinus heart 
rhythm with premature electrical stimuli, resulting in ventricular fibrillation in the test 
system. Treatment candidates are administered to the test system undergoing 
ventricular fibrillation and treatment candidates which modulate ventricular 
fibrillation are identified as therapeutic strategies for treatment of ventricular 
fibrillation. 

[0012] The present invention also relates to a method of evaluating the 

predisposition of a subject for the induction of ventricular fibrillation from a condition 
of ventricular tachycardia. This method involves providing a subject in ventricular 
tachycardia and monitoring the electrical stimuli in the heart of the subject. It is then 
determined if a sequence of rest interval values correspond to rest interval values 
predicted to lead to ventricular fibrillation consistent with the histogram of Figure 7. 
[0013] Another aspect of the present invention relates to a method of 

identifying treatment candidates as therapeutic strategies for preventing ventricular 
tachycardia from developing into ventricular fibrillation. This involves providing a 
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subject in ventricular tachycardia and monitoring the electrical stimuli in the heart of 
the subject. It is then determined if the initial 3 stimuli in groups of 4 stimuli 
correspond to rest interval values predicted to lead to ventricular fibrillation consistent 
with the histogram of Figure 7. Treatment candidates which prevent occurrence of a 
fourth stimuli corresponding to a rest interval value predicted to lead to ventricular 
fibrillation consistent with the histogram of Figure 7 are identified. Treatment 
candidates which prevent ventricular tachycardia from becoming ventricular 
fibrillation are identified as therapeutic strategies for preventing ventricular 
fibrillation. 

[0014] One of the hypotheses to account for the apparent link between spiral 

waves and VF is the restitution hypothesis, which proposes that VF is caused by the 
breakup of a single spiral wa\Q into multiple self-perpetuating wavelets (Weiss et aL, 
''Chaos and the Transition to Ventricular Fibrillation: a New Approach to 
Antiarrhythmic Drug Evaluation," Circulation 99:2819-2826 (1999); Gilmour et al., 
"Electrical Restitution, Critical Mass, and the Riddle of Fibrillation/' J Cardiovasc 
Electrophysiol 10:1087-1089 (1999); Garfinkel et aL, "Preventing Ventricular 
Fibrillation by Flattening Cardiac Restitution," Proc Natl Acad Sci USA 97:6061- 
6066 (2000); Riccio et al., "Electrical Restitution and Spatiotemporal Organization 
During Ventricular Fibrillation," Circ Res 84:955-963 (1999); Panfilov et aL, 
"Ventricular Fibrillation: Evolution of the Multiple Wave Hypothesis," Philos Trans 
R Soc LondSerA 359:1315-1325 (2001); Fenton et aL, "Multiple Mechanisms of 
Spiral Wave Breakup in a Model of Cardiac Electrical Activity," Chaos 12:852-892 
(2002); Karma, A., "Spiral Breakup in Model Equations of Action Potential 
Propagation in Cardiac Tissue,'' Physc Rev Lett 71:1 103-1106 (1993); and Karma, A,, 
"Electrical Altemans and Spiral Wave Breakup in Cardiac Tissue," Chaos 4:461-472 
(1994), which are hereby incorporated by reference in their entirety). This 
mechanism is similar to that proposed decades ago by Moe (Moe, G., "A Computer 
Model of Atrial Fibrillation," //ear/ J 67:200-220 (1964), which is hereby 
incorporated by reference in its entirety), but with one important difference: it may 
occur in intrinsically homogeneous cardiac tissue. The transition from normal planar 
wave excitation to a single spiral wave and ultimately to multiple wavelets is thought 
to underlie the transition from normal sinus rhythm to ventricular tachycardia and VF 
characteristic of patients who succumb to sudden death (Figure 1). 
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[0015] The exact mechanism for the breakup of spiral waves is unknown. 

However, there is considerable evidence that breakup is closely related to APD 
restitution, wWch is the relationship between action potential duration (APD) and 
diastolic interval (DI, the time interval between action potentials), where APD is 
determined by the preceding DI (Figure 2). During pacing of cardiac tissue at 
progressively shorter cycle lengths, APD decreases, until at sufficiently short cycle 
lengths a period doubling bifurcation occurs and APD begins to alternate between a 
long duration and a short duration, a phenomenon known as APD altemans (Figure 2) 
(Moe, G., **A Computer Model of Atrial Fibrillation," Heart J 67 aOO-llO (1964) 
(Guevara et al., "Electrical Altemans and Period Doubling Bifurcations," IEEE Camp 
Cardiol 562:167-170 (1984); Chialvo et al., "Low Dimensional Chaos in Cardiac 
Tissue," Nature 343:653-657 (1990); Karagueuzian et al., "Action Potential Altemans 
and In'egular Dynamics in Quinidine-Intoxicated Ventricular Muscle Cells. 
Implications for Ventricular Proarrhythmia," Circulation 87:1661-1672 (1993); 
Nolasco et al., "A Graphic Method for the Study of Alternation in Cardiac Action 
Potentials," JAppl Physiol 25:191-196 (1968); Rosenbaum et aL, "Electrical 
Altemans and Vulnerability to Ventricular Arrhythmias," New Engl J Med 330:235- 
241 (1994); Pastore et al., "Mechanism Linking T-Wave Altemans to the Genesis of 
Cardiac Fibrillation," Cir Res 99:1385-1394 (1999), which are hereby incorporated by 
reference in their entirety). The development of APD altemans requires that the slope 
of the restitution relation exceed 1 , whereas if the slope of the restitution is less than 
1 , altemans will not occur. 

100161 The slope of the restitution relation, and the corresponding presence or 

absence of APD altemans, has been linked to spiral wave stability. If the slope of the 
APD restitution relation is < 1, a spiral wave tends to stabilize and produce a periodic 
rhythm, the manifestation of which might be monomorphic ventricular tachycardia. 
If, however, the slope of the APD restitution relation is > 1, a single spiral wave may 
disintegrate into many spiral waves, manifest as VF. These transitions are illustrated 
in Figure 3, using results generated by a computer model. Once initiated, a single 
stable spiral wave can be destabilized by increasing the slope of the restitution 
relation from < 1 to > 1 . Similarly, multiple wavelets can be induced to coalesce into 
a single spiral wave by reducing the slope of the restitution relation from > 1 to < 1 . 
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Thus, a reduction of the slope of the restitution relation from > 1 to < 1 is expected to 
prevent the induction of VF and to convert existing VF into a periodic rhythm. 
[0017] Several recent studies have provided experimental and theoretical 

support for a causal relationship between APD restitution and VF (Garfmkel et ah, 
"Preventing Ventricular Fibrillation by Flattening Cardiac Restitution," Proc Natl 
AcadSci USA 97:6061-6066 (2000); Riccio et al., "Electrical Restitution and 
Spatiotemporal Organization During Ventricular Fibrillation," Circ Res 84:955-963 

(1999) ; Koller et aL, "Dynamic Restitution of Action Potential Duration During 
Electrical Altemans and Ventricular Fibrillation,'Mm JP/iy^y/o/ 275:H1635-H1642 
(1998); Koller et al., ''Effects of [K+]o on Electrical Restitution and Activation 
Dynamics During Ventricular Fibrillation/' Am J Physiol Heart Circ Physiol 
279:H2665-H2672 (2000); Fox et al., ''Spatiotemporal Transition to Conduction 
Block in Canine Ventricle," Circ Res. 90:289-296 (2002); Chen et aL, "Spirals, 
Chaos, and New Mechanisms of Wave Propagation," Pacing Clin Electrophysiol 
(PAC) 20:414-421 (1997); and Panfilov, A,, "Spiral Breakup as a Model of 
Ventricular Fibrillation," Chaos 8:57-64 (1998), which are hereby incorporated by 
reference in their entirety. If the APD restitution relation contains a region of slope > 
1 , APD altemans and VF are induced by pacing at short cycle lengths. If the slope of 
the APD restitution relation is reduced to less than one, APD altemans is suppressed 
and VF is not induced. Furthermore, if a fibrillating ventricle is exposed to a drug 
that reduces the slope of the restitution to less than 1 , VF is converted to a periodic 
rhythm sustained by a single stable spiral wave (Garfinkel et al., "Preventing 
Ventricular Fibrillation by Flattening Cardiac Restitution," Proc Natl Acad Sci USA 
97:6061-6066 (2000); Riccio et al., "Electrical Restitution and Spatiotemporal 
Organization During Ventricular Fibrillation," Circ Res 84:955-963 (1999); and 
Koller et al., "Effects of [K+]o on Electrical Restitution and Activation Dynamics 
During Ventricular Fibrillation," Am J Physiol Heart Circ Physiol 279:H2665-H2672 

(2000) , which are hereby incorporated by reference in their entirety). 
[0018] Unfortunately, interventions that suppress VF experimentally 
(verapamil (Riccio et al., "Electrical Restitution and Spatiotemporal Organization 
During Ventricular Fibrillation," Circ Res 84:955-963 (1999), which is hereby 
incorporated by reference in its entirety), bretylium (Garfinkel et aL, "Preventing 
Ventricular Fibrillation by Flattening Cardiac Restitution," Proc Natl Acad Sci USA 
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97:6061-6066 (2000), which is hereby incorporated by reference in its entirety) and 
hyperkalemia (Keller et al., "Effects of [K+]o on Electrical Restitution and Activation 
Dynamics During Ventricular Fibrillation/' Am J Physiol Heart Circ Physiol 
279:H2665-H2672 (2000), which is hereby incorporated by reference in its entirety) 
have unwanted actions that severely limit their clinical utility. Nevertheless, the 
effects of these interventions on restitution and VF have provided valuable insights 
regarding new, potentially more clinically relevant, drug targets. In particular, 
attempts to understand the mechanism by which calcium channel blockers alter APD 
restitution have led to the realization that increasing selected outward repolarizing 
currents also may flatten restitution, as presented in the examples. 
[0019] The physiological mechanism associated with the initiation of 

ventricular fibrillation are broadly debated. The present invention affords the ability 
to characterize and identify the events initiating ventricular fibrillation. As 
demonstrated by the examples, the present invention also affords novel methods of 
identifying candidates for the treatment and prevention of ventricular fibrillation. 



BRIEF DESCRIPTION OF THE DRAWINGS 

[0020] Figure 1 shows the proposed evolution of cardiac wave propagation 

patterns underlying the transition from a normal cardiac rhythm to ventricular 
fibrillation. 1 ransition from normal sinus rhythm (NSR) to ventricular tachycardia 
(VT) and ventricular fibrillation (VF), as recorded on the surface electrocardiogram 
are shown. Corresponding transitions from a planar wave to a single wave to multiple 
spiral waves are shown. 

[0021] Figures 2A-C show the restitution of action potential duration (APD). 

Figure 2 A shows the definition of APD, diastolic interval (DI), and cycle length (CL). 
The duration of the n+lth action potential (APDn+i) is a function of the preceding DI 
(Din). Figure 2B shows the plot of APD versus preceding DI. Figure 2C shows the 
plot of APD versus CL. Crosses are APDs during the CLs in which no alternans 
occurs. During APD alternans, unfilled circles are APDs during the long action 
potential, and filled circles are APDs during the short action potential. Modified, with 
permission, from Koller et al. (Koller et al., "'Dynamic Restitution of Action Potential 
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Duration During Electrical Altemans and Ventricular Fibrillation," ^/w J Physiol 
275:H1635-H1642 (1998), which is hereby incorporated by reference in its entirety). 
[0022] Figures 3 A-B show the wave patterns in a computer model of cardiac 

activation. Sequential snapshots (A-1) of activation during the development of a 
stable single spiral wave under conditions in which the slope of the restitution relation 
for action potential duration is < 1 (Figure 3 A). The transition form a single spiral 
wave into many spiral waves for a restitution slope >1 (A-O, with some intervening 
steps deleted: Figure 3B). 

[0023] Figure 4A shows a schematic representation of the relationships 

between stimulus period (T), action potential duration (D), diastolic interval (I) and 
memory (M). Figure 4B shows an example of the stimulation protocol. The model 
fibre was paced for ten beats at a constant SI SI S'l-S'l interval, after which a series 
of premature stimuli ( S2 - 85 S'2-S'5 ) was delivered. R indicates the refractory period 
of a cardiac cell, during which the delivery of a stimulus does not produce a 
propagated response. 

[0024] Figure 5 shows the spatial dispersion of repolarization induced in an 

ionic model of a homogenous 1 -dimensional cable by multiple premature stimuli. 
Note that short duration action potentials at the site of stimulation increase in duration 
as they propagate alone the fiber, whereas long duration action potentials shorten. 
The resultant discordant altemans pattern culminates in conduction block. 
[0025] Figure 6 shows an example of a type II conduction block in the 

coupled maps model. Each panel shows I as a function of cell number for the last S; 
[[S]]- 500 ms beat (left panel) and for beats S:2 [[S2]]- 207 ms, SH [[S3]]- 1 46 ms, 
SM [[S4]J- 148 ms and SIS [[S5]]- 106 ms. Note the difference in scale for each 
beat. Conduction block, corresponding to I < I min , occurs near cell 100 on the S'5 
[[S5J] beat and is denoted by two horizontal black lines. 

[0026] Figure 7 shows a histogram of the incidence of type II conduction 

block in the wild type model, a = 0.2, xo ^ 28 and p = 1 7.408. The height at each 
point in (S'2 [[S2]], S13 [[S3]]) space corresponds to the number of blocks found at 
beat SIl [[S3]], S14 [[S4]] or 815 [[S5]] for a given (812 [[S2]], 813 [[S3]]) pair. 
Total counts = 67 205. S'2rmn [[S2min ]]and S'3mm [[S3min ]] are the minimum time 
intervals at points S^ [[S2]] and S13 [[S3]], respectively, that allow conduction down 
the fiber, td is the time interval between activations used to adjust dependence of the 
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D recovery function on 1. P is the velocity variable which was varied to adjust the 
value of V at I = Imin- a is the memory variable (which varied between zero and unity) 
which was used to determine the influence of memory on D. 

[0027] Figures 8A-C show the role of V restitution in type II block. Figure 

8A shows V restitution for four cases: p = 17,408 (y-axis - second line from top; wild 
type), p = 7.704 (y-axis - third line from top; V(Imm) = 50% of Vmax), P 2.2028 (y- 
axis - fourth line from top; V(Imin) ^ 25% of Vmax), and P == 30.236 (y-axis - first line 
from top; V(Imin) = 90% of Vmax)- Figure 8B shows a histogram for V(Imin) = 50% of 
Vmax. Total counts = 8879 (cf 67 205 for wild type). Figure 8C shows a histogram 
for V(Imin) 90% of Vmax- Total counts — 86 032. There were no counts tor V(Imin) ~~ 
25% of Vmax. 

[0028] Figures 9A-D show the role of D restitution in type II block. Figure 

9A shows steady-state D restitution for four cases: xo 28 (y-axis - third line from 
top; wild type; max, slope = 1.09), td = 32 (y-axis - second line from top; max. slope 
= 0.97), td= 40 (y-axis - first line from top; max. slope = 0.82) and td== 20 (y-axis - 
fourth line from top; max. slope = 1.45). Figure 9B shows a histogram for Td = 32. 
Total counts = 39 688. Figure 9C shows a histogram for td 40. Total counts = 21 
185. Figure 9D shows a histogram for td = 20. Total counts = 364 301 . 
[0029] Figures 1 OA-D show the role of M in type II block. Figure 1 OA shows 

the steady-state D restitution for four cases: a = 0.2 (y-axis - second line firom top; 
wild type; max. slope = 1.09), a = 0.58 (y-axis - third line from top; max. slope = 
1 .62), a ^ .08 (y-axis - fourth line from top; max. slope = 2.65) and a = 0. 1 (y-axis - 
first line from top; max. slope = 1 .09). Figure lOB shows a histogram for a = 0.58. 
Total counts 48 387. Figure 1 OC shows a histogram for a ^ 0.8. Total counts 43 
715. Figure lOD shows a histogram for a = 0.1. Total counts =^ 87 433, 
[0030] Figure 1 1 shows examples of MAP recordings obtained from the right 

(RV) and left (LV) ventricles of a closed-chest anesthetized beagle during delivery of 
multiple premature stimuli (CLvf) to the RV. CLvf were chosen to produce maximum 
dynamically-induced dispersion of repolarization, as predicted from the results of the 
coupled maps model shown in Figure 7. 

[0031] Figure 12 shows the incidence of conduction block in Beagles. The 

incidence of conduction block followed the delivery of 4 premature stimuli in the 
coupled maps memory model (Figure 7), and following the delivery of CLvf in 



12 



anesthetized beagles (filled circles). CLvpinduced VF in all 8 dogs. In contrast, 
patterns of cycle lengths that fell outside CLvf did not induce VF (unfilled circles). 
[0032] Figure 1 3 shows the spontaneous non-sustained polymorphic 

ventricular tachycardia in a German Shepherd dog. 

[0033] Figures 14A-C show the effects of verapamil on action potential 

duration (APD) restitution and ventricular fibrillation (VF). Figure 14A shows an 
example of the reduction in the slope of the APD restitution relation after exposure to 
verapamil. Figure 14B shows the corresponding reduction in the magnitude of APD 
altemans (AM) during pacing at different basic cycle lengths (BCL) after exposure to 
verapamil. Figure 14C shows the effects of verapamil on the composite Fast Fourier 
Transform (FFT) of monophasic action potential (MAP) recordings during VF in 
arterially perfused canine ventricle. Verapamil, added after 35 min of VF had elapsed, 
converted the broad FFT spectra during VF to a single frequency. Modified, with 
permission, from Riccio et al, (Riccio et al., "Electrical Restitution and 
Spatiotemporal Organization During Ventricular Fibrillation," Circ Res 84:955-963 
(1999), which is hereby incorporated by reference in its entirety). 
[0034] Figures 1 5 A-B show the putative mechanism for APD altemans and 

the suppression of alternans by calcium channel blockade. Under control conditions, 
alternation of the length of the diastolic intervals (Dl) during APD alternans causes 
alternation of the magnitude of Ica, secondary to incomplete recovery of Ica from 
inactivation during the short DI (Figure 1 5 A), Reduction of Ica decreases APD, 
which prolongs DI, so that DI following each APD is sufficiently long to allow 
complete recovery of Ica (Figure 15B). In the absence of alternations of Ica, APD 
altemans in suppressed. 

[0035] Figures 16A-D show the effects of HERG overexpression on iRr and 

APD alternans in isolated canine cardiac myocytes. Current traces of Ikf (Figure 16A) 
and action potentials during pacing at CL ^ 200 msec (Figure 16B) were obtained 
from myocytes in the day of isolation (Day 0). Current traces of Ikf (Figure 16C) and 
action potentials during pacing at CL = 200 msec (Figure 16D) were obtained from 
myocytes after 3 days in cell culture and infection with adenovirus (Day 3 + HERG). 
HERG overexpression increased iKr and suppressed APD altemans. 
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DETAILED DESCRIPTION OF THE INVENTION 
[0036J The present invention relates to a method of evaluating the effect of a 

physiological condition on the occurrence of ventricular fibrillation. This method 
involves providing a test system and initiating a ventricular fibrillation inducing 
sequence in the test system by interrupting normal sinus heart rhythm with premature 
electrical stimuli. Following the initiation of ventricular fibrillation a ventricular 
recovery sequence is initiated. The test system is subjected to a physiological 
condition before, during, and after initiating the ventricular fibrillation recovery 
sequence and the physiological conditions which affect ventricular fibrillation 
recovery are identified. 

[0037] In carrying out this aspect of the present invention, the physiological 

condition to be tested can be velocity restitution, action potential duration restitution, 
or cardiac memory. 

[0038] The step of initiating a ventricular fibrillation inducing sequence is 

carried out by: (1) initiating 4 prematvire stimuli under conditions effective to initiate 
an excitatory wave; (2) altering velocity recovery function to slow conduction at 
shorter rest intervals; (3) altering action potential duration recovery function to reduce 
sensitivity of action potential duration to a preceding rest interval; or (4) is carried out 
by increasing cardiac memory. The test system can be a test animal, a tissue, a cell 
culture, or an in vitro system. 

[0039] Another aspect of the present invention relates to a method of 

identifying treatment candidates as therapeutic strategies for the prevention of 
ventricular fibrillation. This method involves providing a test system and initiating a 
ventricular fibrillation inducing sequence in the test system by interrupting normal 
sinus heart rhythm with premature electrical stimuli before, during, and after 
administering the treatment candidate to the test system. Treatment candidates wWch 
prevent the initiation step Irom inducing ventricular fibrillation are identified as 
therapeutic strategies for prevention of ventricular fibrillation. 
[0040] In carrying out this aspect of the present invention, the treatment 

candidate is a pharmaceutical compound, such as a calcium channel antagonist. 
Alternatively, the treatment candidate is one or more electrical impulses. 
[0041] The step of identifying identifies treatment candidates which achieve 

velocity restitution values consistent with the histograms of Figures 8B and 8C as 
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therapeutic strategies for prevention of ventricular fibrillation. Alternatively, the step 
of identifying identifies treatment candidates which achieve potential duration 
restitution values consistent with the histograms of Figures 8B and 8C as therapeutic 
strategies for prevention of ventricular fibrillation. Finally, the step of identifying 
identifies treatment candidates which achieve cardiac memory values consistent with 
the histograms of Figures 9B and 9C as therapeutic strategies for prevention of 
ventricular fibrillation. 

[0042] Again, the test system can be a subject (i.e., a test animal or patient), a 

tissue, a cell culture, an in vitro system. The step of initiating a ventricular fibrillation 
inducing sequence can be carried out by initiating 4 premature stimuli under 

conditions effective to initiate an excitatory wave. 

[0043] A further aspect of the present invention relates to a method of 

identifying treatment candidates as therapeutic strategies for treating ventricular 
fibrillation. This method includes providing a test system and initiating a ventricular 
fibrillation inducing sequence in the test system by interrupting normal sinus heart 
rhythm with premature electrical stimuli, resulting in ventricular fibrillation in the test 
system. Treatment candidates are administered to the test system undergoing 
ventricular fibrillation and treatment candidates which modulate ventricular 
fibrillation are identified as therapeutic strategies for treatment of ventricular 
fibrillation. 

[0044] In carrying out this aspect of the present invention, the treatment 

candidate is a pharmaceutical compound, such as a calcium channel antagonist. 
Alternatively, the treatment candidate can be one or more electrical impulses. 
[0045] The step of identifying identifies treatment candidates which achieve 

velocity restitution values consistent with the histograms of Figures 8B and 8C as 
therapeutic strategies for prevention of ventricular fibrillation. Alternatively, the step 
of identifying identifies treatment candidates which achieve potential duration 
restitution values consistent with the histograms of Figxores 9B and 9C as therapeutic 
strategies for prevention of ventricular fibrillation. As yet another altemative, the step 
of identifying identifies treatment candidates which achieve cardiac memory values 
consistent with the histograms of Figures lOB and IOC as therapeutic strategies for 
prevention of ventricular fibrillation. 
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[0046] In carrying out this aspect of the present invention, the test system can 

be a subject (i.e. a test animal or patient), a tissue, a cell culture, or an in vitro system. 
The step of initiating a ventricular fibrillation inducing sequence is carried out by 
initiating 4 premature stimuli under conditions effective to initiate an excitator>' wave, 
[0047] The present invention also relates to a method of evaluating the 

predisposition of a subject for the induction of ventricular fibrillation from a condition 
of ventricular tachycardia. This method involves providing a subject in ventricular 
tachycardia and monitoring the electrical stimuli in the heart of the test animal. It is 
then determined if a sequence of rest interval values correspond to rest interval values 
predicted to lead to ventricular fibrillation consistent with the histogram of Figure 7. 
[0048] Another aspect of the present invention relates to a method of 

identifying treatment candidates as therapeutic strategies for preventing ventricular 
tachycardia from developing into ventricular fibrillation. This involves providing a 
subject in ventricular tachycardia and monitoring the electrical stimuli in the heart of 
the test animal. It is then determined if the initial 3 stimuli in groups of 4 stimuli 
correspond to rest interval values predicted to lead to ventricular fibrillation consistent 
with the histogram of Figure 7. Treatment candidates which prevent occurrence of a 
fourth stimuli corresponding to a rest interval value predicted to lead to ventricular 
fibrillation consistent with the histogram of Figure 7 are identified. Treatment 
candidates which prevent ventricular tachycardia from becoming ventricular 
fibrillation are identified as therapeutic strategies for preventing ventricular 
fibrillation. 

|0049] In carrying out this aspect of the present invention, the treatment 

candidate is a pharmaceutical compound, such as a calcium channel antagonist. 
Alternatively, the treatment candidate can be one or more electrical impulses. 
[0050] In one embodiment of this aspect of the present invention, the 

identifying step identifies treatment candidates which achieve velocity restitution 
values consistent with the histograms of Figures 8B and 8C as therapeutic strategies 
for prevention of ventricular fibrillation. Alternatively, the step of identifying 
identifies treatment candidates which achieve potential duration restitution values 
consistent with the histograms of Figures 9B and 9C as therapeutic strategies for 
prevention of ventricular fibrillation. A further alternative is that the identifying step 
identifies treatment candidates which achieve cardiac memory values consistent with 



16 



the histograms of Figures lOB and IOC as therapeutic strategies for prevention of 
ventricular fibrillation. 

[0051] Previous theoretical and experimental studies of the impact of dynamic 

heterogeneity of repolarization on the development of conduction block were 
conducted in computer models of a homogenous 1 -dimensional cable and in isolated 
Purkinje fibers during fixed pacing. As the pacing cycle length was shortened 
progressively, a transition from 1:1 phase locking to a concordant altemans (2:2 
locking) to discordant altemans to intermittent 2:1 conduction block occurred (Fox et 
al, ''Spatiotemporal Transition to Conduction Block in Canine Ventricle," Circ Res, 
90:289-96 (2002), which is hereby incorporated by reference in its entirety). These 
results confirmed the predictions of previous computer modeling studies of these 
phenomena (Watanabe et al., "Mechanisms for Discordant Altemans," J Cardiovasc 
Electrophysiol 12:196-206 (2001), which is hereby incorporated by reference in its 
entirety) and extended experimental observations that the development of discordant 
altemans was associated with local conduction block and the initiation of reentry 
(Pastore et al., "Mechanism Linking T-Wave Altemans to the Genesis of Cardiac 
Fibrillation," Circulation 99:1385-94 (1999); Pastore et al., "Role of Structural 
Barriers in the Mechanism of Altemans-Induced Reentry," Circ Res. 87:1 157-63 
(2000), which are hereby incorporated by reference in their entirety). In addition, 
they demonstrated that complex dynamics was possible in a homogenous 1 - 
dimensional cable, in the absence of intrinsic heterogeneity or anisotropy, and that the 
mechanism for the dynamics could be described using simple unidimensional maps 
for the restitution of APD and CV. 

[0052 1 The evolution of dynamics during fixed pacing occurs relatively 

gradually, which allows this phenomenon to be studied in the quasi-steady-state. As a 
result, a detailed description of the dynamics and identification of the underlying 
mechanism is possible. It is unlikely, however, that prolonged activation at a short 
constant cycle length initiates VF in patients at risk for sudden death. Therefore, 
investigations into a scenario that is more common in such patients - initiation of VF 
following a run of non-sustained VT, were begun. In the present invention, computer 
models of a homogenous 1 -dimensional cable were generated using a coupled maps 
memory model (Fox et al., "Period-Doubling Instability and Memory in Cardiac 
Tissue," Phys Rev Lett 89:138101 (2002), which is hereby incorporated by reference 



17 



in its entirety) and an ionic model (Fox et al., "Ionic Mechanism of Electrical 
Altemans," Am J Physiol Heart Circ Physiol 282:H5 16-30 (2002), which is hereby 
incorporated by reference in its entirety). The models were paced at a cycle length of 
1000 or 500 ms for 50 beats, after which a series of premature stimuli was delivered 
(Figures 4 and 5). Premature stimuli delivered at short DI (e.g., SI2 [[S2]]in Figure 5) 
elicited action potentials at the site of stimulation that, in accordance with the relevant 
APD and CV restitution relations were short in duration and conducted slowly. Slow 
conduction along the fiber provided progressively more recovery time as the action 
potential conducted down the fiber. Accordingly, DI increased along the fiber, as did 
the corresponding APD, resulting in a discordance between APD at the site of 
stimulation and at the distal end of the fiber. Conversely, premature stimuli delivered 
at longer DI (S13 [[S3]] in Figure 5) produced longer duration action potentials that 
conducted more rapidly. The more rapidly conducting action potentials impinged 
progressively on the preceding DI as they conducted along the fiber, creating a shorter 
DI and corresponding APD. Once again, discordant APD between the site of 
stimulation and the distal end of the fiber resulted. The spatial dispersion of 
repolarization caused by this sequence of cycle lengths could be amplified by the 
delivery of multiple premature stimuli, eventually resulting in conduction block ( S'5 
[[SslJin Figures). 

[0053] To more completely characterize this sequence of events, all possible 

combinations of SVS'^ [[Sa-Ssjjthat produced conduction block in the coupled maps 
model (block also occurred with many fewer combinations of S'7"S% [[S7-S^]] and 
many more combinations of S'^^S' ? [[S2-S7]]) were determined. Although 
approximately 67,000 combinations produced conduction block, the distribution of 
such combinations was not random, but was clustered according to specific S^ [[S2]] 
and SI3 [[S3]] interv^als, as shown in Figure 7. The incidence of block was greatest 
when S^ [[S2]] was delivered at a coupling interval near S^iS^min [[SiS2mjn]] (the 
shortest SlxSl2_[[SiS2]] interval that produced a propagated response) and S^ [[S3]] 
was delivered at a coupling interval approximately 50 ms longer than S'^SNmin 
[[SaSBmin]]' Although not shown on this 3-dimensionaI plot, the SI3SI4 [[S3S4]] 
interval most often associated with block was near S^Sl4min [[S3S4min]]? whereas the 
SliSls [[S4S5]] interval was approximately 10 ms longer than S^4S^5min [[S4S5min]]' 
This short-long-short-intermediate pattern of stimulation produced discordant 
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altemans and conduction produced discordant altemans and conduction block of the 
type shown in Figure 5, Also determined was the impact of changes in the slope of the 
APD restitution relation, the magnitude of memory, the shape of the CV restitution 
relation and the strength of cell coupling. In general, the incidence of conduction 
block was reduced by decreasing the slope of the APD restitution relation, increasing 
memory, increasing cell coupling and extending the CV restitution relation to permit 
slow conduction (thereby preventing "head-tail" interactions) (Fox et al.. Dynamic 
Mechanism for Conduction Block in Heart Tissue," New J Phys, 5:101.1-101.14 
(2003), which is hereby incorporated by reference in its entirety). 

EXAMPLES 

[0054] The following examples are provided to illustrate embodiments of the 

present invention but are by no means intended to limit its scope. 

Example 1 - Methods of Examining the Dynamic Mechanism for Conduction 
Block in Heart Tissue 

[0055] The mechanism for conduction block following the delivery of 

multiple premature stimuli using a one-dimensional computer model of a canine heart 
fiber is presented. The fiber was stimulated at one end at a cycle length near normal 
canine sinus rhythm (SJ, [[SI]] 500 ms). Atlter ten beats at S:i [[SI]], four 
premature stimuli were delivered (S:2[[S2]], S:3[[S3]], S:4[[S4]] and S:5[[S5]]). The 
coupling intervals between these premature intervals were varied and conduction of 
the resultant action potentials was observed (Figure 4). For each combination of 
premature stimuli, one of three possible outcomes occurred: 1) the stimulus elicited an 
action potential that propagated down the entire fiber; 2) the stimulus did not produce 
an action potential (Type I block); 3) the stimulus elicited an action potential that 
blocked before reaching the end of the fiber (Type II block). 
[0056] The first two cases would not be conducive to the development of 

wave break and initiation of reentry, in that conduction of the premature response 
either does not occur at all or occurs equally well everywhere along the fiber. The 
third case, however, could lead to wave break and spiral wave initiation if it occurred 
in two- and three-dimensional tissue, provided the block was local. The development 
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of local block would be facilitated in intact myocardium by twist anisotropy and 
intrinsic heterogeneity, 

[0057] To assess the vulnerability of the simulated tissue to Type II block, the 

S'3-SM-S'5 f[S2-S3-S4-S5n combinations were applied in the following way: the 
S'1"S'2 [[S1-S2]] interval was varied from the minimum value that conducted (S'l- 
S:2xnm_[[Sl-S2min]]) to 1 ■S^2min IfS 1 ■S2n.inl1 + 20 ms. For each Sl^S!! [[S1-S2]] 
interval, the S^2^S^3 [[S2-S3]] interval and the S'3-S'4 [[S3-S4]] inter\^al were varied 
in combination from the minimum value that conducted up to a value of 250 ms for 
each interval. Finally, for each S^2-S^3-S'4 [[S2-S3-S4]] combination that 
conducted, the S'4-S'5 [[S4-S5]] interval was varied from the minimum that 
generated an action potential at the site of stimulation to the minimum value that 
conducted down the entire fiber. If no S^2-S'3-SM>S^5 |TS2^S3-S4-S511 combination 
was found that produced conduction block for more than a 20 ms window in either the 
S^2-S^3 [[S2-.S3]] interval or the S^3-SM [[S3-S4]] interval, the search was halted for 
that interval. All intervals were varied in steps of 1 ms. 

Example 2 - Coupled Maps Model 

[0058] One aspect of the present invention is conducted using a coupled maps 

model of a one-dimensional cardiac fiber that has been described in detail elsewhere ( 
Fox et al., "Period-Doubling Instability and Memory in Cardiac Tissue," Phys Rev 
Lett 89:138101 (2002) and Fox et al, "Spatiotemporal Transition to Conduction Block 
in Canine Ventricle," Circ Res. 90:289-96 (2002), which are hereby incorporated by 
reference in their entirety). Briefly, the model is based on the equation 

/„,,(x,)=r„,,(;r,)-A,.,(^,). (1) 



Tn+/ (Xi) was the time interval between activations of site It was determined by 
including the time delays caused by the propagation from the pacing site to site , 
which yielded 
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T was the time interval between activations applied to the pacing site and Ax ^ O.I was 
the length of a single cell (time units in ms and space units in mm). The conduction 
velocity Vn(xi) depended only on / through the velocity recovery function K„ = c(I^ 
given by c(I)'==Vmax(l-exp(-(I+fi)/S)). V^ax^^^O. 72, 6=^14, and /? was varied to adjust the 
value of F at l=Imm^ D was determined locally based on a memory model mapping (8) 
given by: 



Tm was the time constant of accumulation and dissipation of memory (both constants 
were chosen to be the same). A = 88, B = 722, C = 40^ and = 180 (Fox et aL, 
"Period-Doubling Instability and Memory in Cardiac Tissue,'' Phys Rev Lett 
89:138101 (2002), which is hereby incorporated by reference in its entirety), a and 
td were varied in this study, a (which varied between 0 and 1) determined the 
influence of memory on D, and td was used to adjust the dependence of the D 
recovery function on / (Fox et ah, 'Teriod-Doubling Instability and Memory in 
Cardiac Tissue," Phys Rev Lett 89:138101 (2002), which is hereby incorporated by 
reference in its entirety). 

[0059] Coupling between sites was included by using the diffusion terms from 

Echebarria and Karma (Echebarria et al., ''Instability and Spatiotemporal Dynamics of 
Altemans in Paced Cardiac Tissue," Phys Rev Lett 88:208101 (2002), which is hereby 
incorporated by reference in its entirety). These terms modeled the electrotonic 
current that flowed out of (into) a cell and into (out of) its neighbor if the action 
potential of the first cell was longer (shorter) than that of its neighbor. Including 
diffusion then yielded 




(3) 



(4) 
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with 1^=1.0 and w='.35. Discretizing the derivatives produced a tri-diagonal linear 
system of equations that could then be solved easily. The defining equation for the 
model was therefore 



^„.> =/(A^„.,,/„) + #'V^/>„,, -H'V£>„,, (5) 

Ir,Ax.) = r + y — y — D Ax ) 



Finally, conduction block was modeled by setting/ = 0 for / < Imm = 2. 



Example 3 - Characterization of Type II Conduction Block in Wild-type Model 

[0060] Figure 6 shows a typical example of type II conduction block in the 

baseline model (hereafter called the wild type model). Despite the homogeneous state 
that existed at the end of the SH [[SI]] stimuli, dynamical heterogeneity developed 
following delivery of S^2 [[S2]], secondary to F restitution. Heterogeneity was 
manifest as a short-to-long gradient in /, which was magnified by subsequent 
premature stimuli until the action potential generated by the SI5 [[S5]] stimulus 
encountered a region where / < /min and conduction failed. 

[0061] Figure 7 shows a histogram of type II conduction block at various S12 

[[S2]] and SH [[S3]] intervals generated from the wild type model. Note the presence 
of a large peak centered at an S12 [[S2]] very close to the minimum value for 
conduction and at an S13 [[S3]] approximately 50 ms longer than the minimum value 
for conduction. Surrounding the peak is a flat plateau region of five to ten blocks per 
bin that extends throughout the entire (812 [[S2]], SH [[S3]]) region that was explored 
numerically. An overwhelming majority of these instances of block occurred after the 
S'5 [[S5]] stimulus. However, six examples of *early block' were found after an SI4 
[[S4]] stimulus, and one was found after an SH [[S3]] stimulus. The latter explain the 
presence of the plateau in the histogram; the plateau is made up of 'degenerate' counts 
in which long S:2 [[S2]] or SIl [[S3]] intervals that maintained the initial SJ, [[SI]] 
homogeneity were followed by early block intervals leading to block at S15 [[S5]]. 
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[0062] Since conduction block occurs when / < Imm, changes in Imm will 

produce significant changes in the incidence of conduction block. Because Imm is a 
component of the functions that control F restitution, D restitution and M, changes in 
/min influence all three dynamic variables. Imm was fixed and individually the roles that 
F restitution, D restitution and M played in determining the risk of type II block, as 
assessed by the S'2-S^3-SM-S^5 [[82-83-84-85]] protocol, was examined. 

Example 4 ~ Contribution of V Restitution to Conduction Block 

[0063] f'^ restitution played two important roles in the development of 

conduction block. First, as mentioned above, F restitution was the source of the initial 
heterogeneity after delivery of the SI2 [[82]] stimulus. The [[S2]] stimulus was 
delivered just after the refractory period of the last SM [[81]] beat so that / at the site 
of stimulation was short (see Figure 4). The short / produced slow V, which increased 
/ for cells distal to the stimulus site. This initial dispersion in / following §12 [[82]] 
was magnified by steep D restitution at each subsequent beat, ultimately leading to 
conduction block several beats later, 

[0064] 8econd, the value of the T restitution curve at Imm was an important 

parameter in determining the likelihood of conduction block, as illustrated by Figure 
8. Figure 8 A shows F restitution curves with four different values for the V(Imm) (90, 
75, 50, and 25% of Fmax). V(Imm) in the wild type model was 75% of Fmax, Figure SB 
shows a histogram for V(Imm) = 50% of Fmax , and Figure 8C shows the histogram for 
V(Imm) = 90% of Fmax. For V(Imm) = 90% of Fmax, the modcl had roughly 30% more 
instances of conduction block than in the wild type model, whereas the 
model with the lower cut-off values produced dramatically fewer instances of 
conduction block. In fact, the model with V(Imm) = 25% of Fmax produced no cases of 
type II block. This result did not depend on the steepness of the F restitution curve. 
For example, no cases of type II block were found for a model with /3 = -1.4 and 6 = 
2.0, which also had V(Imm) = 25% of Fmax, but had a much steeper slope. The lack of 
type II block when V(Imm) = 25% of Fmax can be understood by noting that if the F 
restitution curve approaches zero at /min, then a propagating wave that approaches a 
region with a very small / can slow down to allow / to increase just in front of it, 
permitting continued conduction. 
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Example 5 — Contribution of D Restitution to Conduction Block 

[0065] As mentioned above, D restitution contributed to conduction block by 

magnifying heterogeneity in /. If the D restitution function had weak dependence on /, 
any initial dispersion of /due to F restitution would be lessened on subsequent 
stimulations. Figure 9 shows the effect of altering D restitution slope by increasing 
and decreasing the parameter td- Figure 9A shows steady state restitution curves for 
four different values of td. The wild type model was td = 28. As shown in Figures 9B 
and 9C, increasing td led to a significantly smaller peak in the histogram, as well as 
fewer total instances of conduction block. Conversely, decreasing td, which led to a 
steeper steady state restitution slope, had significantly more counts (Figure 9D), 
However, it is important to note that even for the td = 40 model, which had a 
maximum steady state slope of 0,82, instances of type II block were found. 

Example 6 - Contribution of M to Conduction Block 

[0066] Because increasing M in the coupled maps model has been shown to 

eliminate steady state altemans, it was hypothesized that increasing A/ also would lead 
to fewer instances of type II block. Figure 10 illustrates the effect of changing Mon 
the likelihood of conduction block in the 101.10 model by increasing and decreasing 
the parameter a. Figure 1 OA shows steady-state restitution for four different values of 
a. The wild type model was a = 0.2. In Figure lOB, a = 0.58, which eliminated 
altemans but maintained a slope ^ 1 . In fact, the steady state restitution slope was 
1 .62, even larger than in the wild type model. In Figure 1 DC, a = 0.8 and the 
maximum steady state restitution slope was 2,65. Both Figures lOB and IOC show 
smaller peaks and fewer total counts than the wild type model. In Figure lOD, M was 
decreased (a = 0. 1) and the maximum steady state slope is 1 ,09. The decreased 
memory model shows more instances of conduction block than the wild type model. 
This result did not depend on the magnitude of D, For example, shifting the curve 
from Figure IOC by adding 70 ms to the function / (Mn+\, In) did not produce more 
instances of type II block (in fact slightly fewer (39 229) were found). Similarly, 
shifting the curve from Figure I0A[[10D]] down by subtracting 70 ms from f(Mn+\, 
In) produced roughly the same number of blocks (94 151) as in Figure lOD. Two 
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other observations can be made. First, increasing Mdid not seem to be as effective as 
increasing td in reducing the likelihood of block. Second, the steady-state D 
restitution function is not predictive of the likelihood of block. 

Example 7 - Contribution of Electronic Interactions to Conduction Block 

[0067] Electrotonic effects have been shown to prevent altemans in some 

models (Cytrynbaum et al., ''Stability Conditions for the Traveling Pulse: Modifying 
the Restitution Hypothesis," Chaos 12:788-99 (2002), which is hereby incorporated 
by reference in its entirety). Therefore, changes in electrotonic current may be 
expected to have an effect on the incidence of conduction block. The coupled maps 
model allowed for the complete removal of electrotonic interactions by eliminating 
the diffusion terms in equation (4). As shown in Figure 6, the potential for conduction 
block was present when the previous beat generated a short-to-long pattern in /. This 
pattern led to a short-to-long pattern in the following D. In the wild type model, both 
the first and second derivative diffusion terms tended to decrease the long D and 
increase the short D, reducing the amount of heterogeneity and thereby reducing the 
likelihood of type II block. Accordingly, elimination of the diffusion terms was 
expected to increase the incidence of type II block and this expectation was 
confirmed; total counts of type II block were 88 603. 

[0068] In the present invention, type II conduction block following the 

delivery of multiple premature stimuli was caused by the same sequence of events 
that caused conduction block during sustained rapid pacing in previous studies (Fox et 
aK "'Spatiotemporal Transition to Conduction Block in Canine Ventricle," Circ Res. 
90:289-96 (2002), which is hereby incorporated by reference in its entirety). 
Stimulation at a short / produced a gradient of increasing / (and, consequently, of D 
and ^0 along the fibre. The succeeding stimulus, if delivered at an appropriate 
interval, encountered a progressively decreasing / as it propagated down the fibre (i.e., 
the new wavefront encountered the waveback of the previous excitation). Thus, a 
short- long sequence of D at the site of stimulation was associated with the 
development of a long— short sequence of D at the opposite end of the fibre, similar to 
the pattern established during stable discordant alternans (Pastore et al., "Mechanism 
Linking T-Wave Altemans to the Genesis of Cardiac Fibrillation," Circulation 
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99:1385-94 (1999); Qu et al., ''Mechanisms of Discordant Altemans and Induction of 
Reentry in Simulated Cardiac Tissue," Circulation 102:1664-70 (2000); and 
Watanabe et al., "Mechanisms for Discordant Alternans," J Car J/ova^x* 
Electrophysiol 12:196-206 (2001), which are hereby incorporated by reference in 
their entirety). When the gradient in / along the fibre was sufficiently steep in the 
ascending direction, the subsequent action potential encountered an / < Imm and 
conduction block occurred. 

[0069] Previous studies have shown that the location of the block can be 

predicted, given the spatial profile of D for the preceding beat (Sampson et al., 
"Simulation and Prediction of Functional Block in the Presence of Structural and 
Ionic Heterogeneity," Am J Physiol 281 :H2597-603 (2001), which is hereby 
incorporated by reference in its entirety). If this information is not known, how^ever, 
then developing a quantitative theory for the type II block observed with present 
invention becomes difficult, not only because the model consists of complicated 
nonlinear difference equations in both space and time, but also because the 
phenomenon of interest is a transient. Therefore, the usual tools of nonlinear 
dynamics, such as finding steady-state solutions and analyzing their stability, do not 
apply. Still, inspection of the defining equations produces a qualitative understanding 
of the results presented in this study. The mechanism for conduction block relies on 
three determinants. First, a non-constant F restitution function produces spatial 
heterogeneity in /. Next, this heterogeneity is amplified (or at least it does not decay 
rapidly) due to sensitive dependence of D on /. Finally, the F restitution function does 
not drop sufficiently close to zero at very short / , causing the wave to collide with the 
refractory period of the preceding excitation. 

[0070] Preventing the first aspect of type II block could only be done by 

forcing the F restitution curve to be constant. Even in that case, type II block could 
still occur if there were intrinsic heterogeneity in the fibre. Therefore, the likelihood 
that type II block could be diminished by two different methods was explored. The 
first method was to alter the F restitution function so that conduction slowed at very 
short /, thus preventing encroachment on the refractory period. It is important to note 
that this method is only viable if the decrease in conduction velocity occurs at short /. 
Decreasing maximum conduction velocity is in fact pro-arrhythmic, because it 
decreases the wavelength (Fenton et al., "Multiple Mechanisms of Spiral Wave 
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Breakup in a Model of Cardiac Electrical Activity," Chaos 12:852-92 (2002) , which 
is hereby incorporated by reference in its entirety). Although selective slowing of 
conduction at short rest intervals may create a greater dispersion of refractoriness after 
a premature beat, the benefits of this intervention seem to outweigh the risks, at least 
under the conditions of this study. Whether a similar strategy will suppress 
conduction block in 2D or 3D tissue with realistic intrinsic heterogeneity remains to 
be determined. 

[0071] The second method for reducing the likelihood of type II block was to 

reduce the sensitivity of D on /. If the system is uni-dimensional, this is done by 
decreasing the D restitution slope. The lower the slope of the D restitution relation, 
the lower the likelihood of spatial dispersion of repolarization secondary to dispersion 
of I, Hov/ever, the model used in the present invention did not have a simple 
relationship between the slope of the steady-state D restitution function and the 
likelihood of block. Instead, the important quantity is the sensitivity of D to the 
preceding / , given by 

(6) 



[0072] The quantity y is the amplification factor relating the dispersion in D to 

the dispersion in the preceding /. If y is small, then any initial dispersion in / will 
decay quickly. If y is large, the dispersion of / can persist or even be amplified. It can 
be shown that at short cycle lengths the first term dominates, giving 



B 



(7) 



[0073] M* and /* are the steady state values at a given cycle length. This 

equation is of limited quantitative use since the system never reaches steady state (and 
the cycle length is continuously changing). Still, it suggests that the likelihood of 
block can be reduced by either increasing To or by increasing a. The simulation results 
support this qualitative argument. This result highlights the fact that the dynamics 
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during premature stimulation cannot be predicted from the steady-state D restitution 
function. In fact, the more relevant experimental measure would be the [fSl- 
S2]] D restitution function at very short 8^2 [[S2]]. The slope of this curve at short / 
would directly measure the dispersion of D after a very short S'l-S'2 [[S1-S2]] 
interval had produced some initial dispersion in /. However, relating the S'l-S'2 
[[S1-S2]] curve to any quantitative prediction about type II block is not 
straightforward . 

[0074j Programmed electrical stimulation has been used extensively in the 

past to induce ventricular arrhythmias in patients thought to be at risk for the 
development of ventricular fibrillation (Anderson et al., "Clinical Value of Cardiac 
Electrophysiology Studies, Cardiac Electrophysiology. From Cell to Bedside 2"^ edn, 
ed D.P, Zipes and J, Jalife (Philadelphia, PA: Saunders) 1 133-50 (1995), which is 
hereby incorporated by reference in its entirety). However, the stimulation protocol 
used in the clinical studies differs significantly from that used here, in that arrhythmia 
induction was attempted using tightly coupled premature stimuli, which produced 
block similar to the type I block. Once the patient's arrhythmia had been induced 
using multiple premature stimuli, the test was repeated after administration of an 
antiarrhythmic drug. In many cases a drug could be found that would suppress the 
patient's inducible arrhythmia. Yet the results of large scale clinical trials eventually 
revealed that patients sent home on such medications were as or more likely to die 
suddenly than patients who did not receive such treatment (Cast, "Preliminary report: 
Effect of Encainide and Flecainide on Mortality in a Randomized Trial of Arrhythmia 
Suppression After a Myocardial Infarction " in New EnglJ Med 321 :406-12 (1989) 
and Waldo et al., "Effect of D-Sotalol on Mortality in Patients with Left Ventricular 
Dysfijnction After Recent and Remote Myocardial Infarction. The SWORD 
Investigators. Survival With Oral D-Sotalol," Lancet 348:7-12 (1996), which are 
hereby incorporated by reference in their entirety). Clearly, the induction of 
arrhythmias using this particular stimulation protocol was not predictive of the 
development of ventricular fibrillation. 

[0075] The failure of the standard method of programmed stimulation to 

accurately assess vulnerability to ventricular fibrillation may relate to the fact that this 
method was designed primarily to induce conduction block and re-entry by perturbing 
intrinsic heterogeneity. The potential contribution of dynamic heterogeneity to the 
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development of conduction block would not be assessed adequately using such a 
protocol. In addition, interventions that were intended to prevent the induction of 
arrhythmias using the standard protocol, such as slowing of conduction or 
prolongation of refractoriness, might be effective in that regard, secondary to a 
reduction in intrinsic heterogeneity, but might not reduce dynamic heterogeneity. 
Consequently, the beneficial effects of such interventions on arrhythmias induced by 
the standard stimulation protocol might not correspond to a reduction in the incidence 
of ventricular fibrillation, if induction of the latter is influenced importantly by 
dynamic heterogeneity. 

[0076] Dynamically induced heterogeneity and conduction block is one of 

many potential mechanisms for spiral break-up in cardiac tissue (see for example 
(Fenton et al., ''Multiple Mechanisms of Spiral Wave Breakup in a Model of Cardiac 
Electrical Activity," Chaos 12:852-92 (2002), which is hereby incorporated by 
reference in its entirety)). Because wave propagation in the heart is extremely 
complex and influenced by many factors, it is unlikely that only one mechanism is 
responsible for induction and maintenance of all cardiac arrhythmias. Nevertheless, 
the results of this study suggest several potential interventions that may reduce the 
probability of arrhythmia induction, including flattening the D restitution slope, 
decreasing V at short / and increasing M There is experimental support for the idea 
that flattening D restitution suppresses ventricular fibrillation (Riccio et aL, 
"Electrical Restitution and spatiotemporal Organization During Ventricular 
Fibrillation," Circ Res 84:955-63 (1 999), which is hereby incorporated by reference in 
its entirety), but altering Kand M with the objective of preventing ventricular 
fibrillation has not been attempted. Tests of these options may not be likely in the near 
future; however, given that drugs that increase the recovery of sodium channels at 
depolarized potentials are not currently available, and alteration of M may not be 
possible until the ionic basis for this phenomenon has been established. Nevertheless, 
judicious alterations of £>, V, and Mmay hold future promise for new and more 
effective therapies for ventricular fibrillation 

Example 8 - Map Recordings of Beagle Right and Left Ventricles 
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[0077] More recently, investigations as to whether dynamically-induced 

spatial dispersion of repolarization promotes the development of conduction block 
and reentrant excitation in settings where underlying intrinsic heterogeneity is present 
were initiated. To that end, the sequence of premature stimulus intervals predicted by 
the computer model to cause conduction block was delivered to the right ventricles of 
closed chest anesthetized normal beagles and MAP recordings were obtained from the 
right and left ventricles. The dogs initially were paced at a cycle length of 400 ms and 
a single premature stimulus (S^2 [[S2]]) was delivered at progressively shorter S'lS'? 
[[S1S2]] intervals until [[S2]] failed to capture (S:iS:2rmn [[SiSsmin]]). The SliSl^ 
[[S1S2]] interval was then set to S^S^rnin [[SiS2min]] + 40 ms, and an SI3 [[S3]] 
stimulus w^as delivered at progressively shorter SI2SI3 [[S2S3]] intervals until Sl 3min 
[[S3mm]] was encountcrcd. The procedure was then repeated to find Sl 4min [[S4min]] 
and Sl smin [[Ssmin]]. Thereafter, a sequence of premature stimuli consisting of 

S:iS:2nun [[SlS2n.in]] + 5 mS, Sl^Slsmin [[S2S3„,in]] + 50 mS, S:3S:4niin [[S3S4min]] + 5 mS 

and Sl4Sl5min [[S4S5min]] + 5-10 ms was delivered. As shown in Figure 1 1, the 
pattems of APD produced by this sequence of premature cycle lengths (CLvf) was 
similar to that produced in the computer models {c.f. Figure 5), Short duration 
responses (S) at the RV recording site were associated with longer duration responses 
(L) at the LV site and vice versa and the pattem of discordant alternans preceded the 
development of VF. 

[0078] The delivery of CLvf has produced VF in all 8 dogs tested thus far, 

whereas the delivery of multiple premature stimuli at sequentially shorter cycle 
lengths has failed to induce VF in any of the dogs (Figure 12). Although, for obvious 
reasons, not all 67,000 possible combinations of cycle lengths predicted to cause 
conduction block in the model, nor more than a few of the intervals that are not 
expected to cause block were tested, the highly reproducible induction of VF by CLvf 
lying near the peak of the plot shown in Figure 1 1 suggests that stimulus protocols 
designed to maximize dynamically-induced spatial heterogeneity of repolarization 
may reliably promote the induction of VF. The exact nature of the spatial dispersion 
produced by such a protocol will require more extensive mapping than can be 
accomplished by two MAP recordings, but the observation that CLvf produces 
discordant altemans between at least 2 RV and LV recording sites suggests that the 
pattern of APD generated by the computer models may pertain to the intact heart. 
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[0079] As a further test of the hypothesis that dynamically-induced spatial 

heterogeneity is arrhythmogenic, it was determined whether a reduction in the slope 
of the APD restitution relation suppresses the development of spatial heterogeneity 
APD and conduction block during premature stimulation, as shown in the computer 
model (Fox et al.. Dynamic Mechanism for Conduction Block in Heart Tissue," New 
J Phys, 5:101.1-101.14 (2003), which is hereby incorporated by reference in its 
entirety). For these studies, the effects of verapamil, which is known to decrease the 
slope of the dynamic restitution relation (Riccio et al., "Electrical Restitution and 
Spatiotemporal Organization During Ventricular Fibrillation," Circ Res 84:955-63 
(1999), which is hereby incorporated by reference in its entirety), in the anesthetized 
closed chest dogs described above (n ^ 8) were tested. MAP duration (MAPD) was 
measured in RV and LV endocardium during programmed stimulation of the RV 
using 4 appropriately timed premature stimuli [[S?-Ss]]). following a train of 

20 stimuli at a constant S^j^S^j^ [[SiSi]] interval of 400 ms. The stimulation protocol 
incorporated the same short-long-short-intermediate sequence shown in Figure 12, 
adjusted as necessary to account for changes in S'2 irn n-S'jmm r[S?mm-S^mmn produced 
by verapamil. VF induction was attempted during control and 30 minutes after i.v. 
verapamil (0.1, 0,3 or 1.0 mg/kg/min), 

[0080] VF induction occurred in all dogs during control and after 0.1 mg/kg 

verapamil and was associated with discordant alternans of MAPD between RV and 
LV. In contrast, after the two higher verapamil doses MAPD alternans was reduced 
in magnitude and was concordant and VF could not be induced. These effects were 
not associated with significant changes in the slope of the standard (S^iSj^ [[8182]]) 
MAPD restitution relation, which was 0.41 ±0,17 during control and was 0.47 + 0.28, 
0.56 + 0.20 and 0.39 ±0.18 after verapamil. Thus, suppression of spatial APD 
heterogeneity by verapamil was associated with an inability to induce VF using 
premature stimuli. However, the doses of verapamil required to produce this effect 
had significant negative dromotropic, chronotropic and inotropic effects: verapamil 
increased the PR interval (fi*om 84.7 ± 3.8 during control to 1 1 1.8 ± 12.2, 132.0 ± 
15.0 and 139.5 ± 18.0 ms after verapamil), decreased heart rate (from 139.0 ± 9.7 to 
1 10.0 ± 2.1, 101.6 ± 3.8, and 95.2 ± 8.8 bpm) and decreased mean blood pressure 
(fi-om 1 16.8 ± 7.3 to 107.4 ± 1 1.4, 108.0 ± 6.3 and 101.2 ± 9.4 mmHg). These results 
suggest that suppression of VT-induced spatial electrical heterogeneity may be one 
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means of preventing VF, but that verapamil is unlikely to be the drug of choice for 
producing such an effect in patients. 

Example 9 - Ventricular Fibrillation in German Shepherd Dogs 

[0081] Although the induction of VF by CLvf in normal beagles confirms the 

predictions made by the computer models, this result is somewhat expected, in that 
the data used to construct the models was obtained from normal canine myocardium. 
This model has the additional limitation of requiring external stimuli to induce VF, 
since these dogs do not have spontaneous ectopy. To test the predictions from the 
computer model more rigorously, next it was investigated whether CLvf induced VF 
in German shepherd dogs that display non-sustained polymorphic VT and have an 
inherited predisposition to sudden death (Figure 13). As discussed above, one of the 
interesting and unanswered questions with respect to these animals is why they often 
live for months with multiple episodes of rapid, polymorphic VT before succumbing 
to sudden death. Based on the response of the beagles to CLvf, it seems possible that 
the German shepherds, despite having copious amounts of ectopy, do not die sooner, 
because they do not routinely generate the proper sequence of coupling intervals, 
insofar as promoting dynamically-induced spatial dispersion of repolarization is 
concerned. 

[0082] As an initial test of this idea, 5 affected German shepherd dogs were 

anesthetized and instrumented as described above for the beagles. CLvf were 
determined from measurements of S^aminiS^5miiu[[S2min-S5min]] and predictions from 
the beagle-based computer model and were delivered to the RV or LV. VF was 
induced in 3 of the dogs, but not in the other two, even after several adjustments of 
the CLvf« Sequences of cycle lengths obtained from Holter recordings during 
multiple episodes of polymorphic VT (CLvt) also were delivered to the dogs. None 
of these sequences induced VF. The sequences during CLvt consisted of 
progressively shorter cycle lengths, in contrast to CLvf, which were short-long-short- 
short. In addition, even during longer runs of non-sustained VT, such as the run 
shown in Figure 13, a CLvf - type pattern did not occur. This observation raises the 
question of why do the German shepherds eventually die suddenly? Unfortunately, 
despite the hundreds of Holter recordings obtained since this colony was established. 
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VF has never been captured. To remedy this situation. Reveal Plus© (Medtronic) 
event recorders were implanted in severely affected dogs. Thus far, the development 
of VF in one dog has been captured. The sequence of cycle lengths preceding VF was 
similar to that used to precipitate VF in the programmed stimulation studies, but 
clearly more episodes of spontaneous VF need to be captured and analyzed before 
definitive conclusions can be drawn. 

Example 10 - Calcium Channel Agonists and Ventricular Fibrillation 

[0083] The suppression of VF by verapamil, in association with a reduction in 

the slope of the APD restitution relation, (Figure 14) suggests that the L-type calcium 
current (Ica) plays a key role in restitution and in the development of VF (Riccio et al., 
''Electrical Restitution and Spatiotemporal Organization During Ventricular 
Fibrillation," Circ Res 84:955-963 (1999) and Chudin et al., "Intracellular Ca^^ 
Dynamics and the Stability of Ventricular Tachycardia," Biophys J 77:2930-2941 
(1 999), which are hereby incorporated by reference in their entirety). To better define 
these roles, an ionic model was developed for the cardiac ventricular myocyte (C VM) 
based on the Winslow (Winslow et aL, "Mechanisms of Altered Excitation- 
Contraction Coupling in Canine Tachycardia-Induced Heart Failure, II.," Circ Res 
84:571-586 (1999), which is hereby incorporated by reference in its entirety) and 
Luo-Rudy models (Luo et al., "A Dynamic Model of the Cardiac Ventricular Action 
Potential, I: Simulation of Ionic Currents and Concentration Changes," Circ Res 
74:1071-1096 (1994), which is hereby incorporated by reference in its entirety) that 
generates physiologically realistic APD altemans over a wide range of pacing cycle 
lengths (Fox et al., "Ionic Mechanism of Cardiac Altemans," Am J Physiol 282:H5 1 6- 
H530 (2002), which is hereby incorporated by reference in its entirety). Using this 
model, the L-type Ca^^ current was implicated as an important determinant of APD 
altemans, according to the mechanism shown in Figure 15. Upon initiation of pacing 
at a short CI-, Ica is fully recovered prior to the first action potential and activates fully 
during the action potential, which results in a long APD. The long action potential 
subsequently is followed by a short diastolic interval, during which Ica fails to recover 
completely from Ca'^'^-induced inactivation. Because of decreased availability of Ica? 



33 



the duration of the subsequent action potential is shorter, which resuhs in a longer 
succeeding diastolic interval, more complete recovery of Ica and a long action 
potential duration. This cycle then repeats, eventually establishing a steady-state 
altemans of Ica and APD. 



Example 11 - Potassium Channel Agonists and Ventricular Fibrillation 

[0084] From the mechanism for APD alternans outlined in Figure 1 5, it was 

anticipated that shortening of APD during rapid pacing will prolong DI, which, in 
turn, might provide adequate time for complete recovery of Ica- If so, reducing the 
magnitude of APD altemans could be accomplished by increasing outward 
repolarizing currents, rather than by decreasing Ica- This idea has been supported by 
additional computer modeling studies in which APD altemans was suppressed by 
increasing any one of several outward repolarizing currents (e.g., the rapid (iKr) and 
slow^ (Iks) components of the delayed rectifier and the inward rectifier (IkO) (Fox et 
al., "Ionic Mechanism of Cardiac Altemans," Am J Physiol 282:H516-H530 (2002), 
which is hereby incorporated by reference in its entirety). Although agonists for these 
currents generally are not available, it might be possible to modify currents such as lFCr 
and Iks by increasing phosphatidyl inositol bisphosphate (PIP2) levels (Bian et al., 
*'HERG K- Channel Activity is Regulated by Changes in Phosphatidyl Inositol 4,5- 
Bisphosphate," Circ Res 89: 11 68-1 176 (2001), which is hereby incorporated by 
reference in its entirety) or by altering the phosphorylation state of the channels 
(Kiehn, J., "Regulation of the Cardiac Repolarizing HERG Potassium Channel by 
Protein Kinase A," Trends Cardiovasc Med 10:205-209 (2001); Heath et aL, 'Trotein 
Kinase C Enhances the Rapidly Activated Delayed Rectifier Potassium Current, l^v. 
Through a Reduction in C-Type Inactivation in guinea Pig Ventricular Myocytes," J 
Physiol 522:391-402 (2000); Marx et al., "Requirement of a Macromolecular 
Signaling Complex for p Adrenegic Receptor Modulation of the KCNQl-KCNEl 
Potassium Channel," Science 295:496-499 (2002), which are hereby incorporated by 
reference in their entirety), provided the channels can be phosphorylated without 
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concomitant phosphorylation of calcium channels or can be upregulated in the 
presence of a calcium channel antagonist, to offset increased Ica secondary to 
phosphorylation (Renter et al., "The Regulation of Calcium Conductance of Cardiac 
Muscle by Adrenaline,' VP/jjy^/o/ 246:49-62 (1977), which is hereby incorporated by 
reference in its entirety). 

[0085] Recently, increasing Ixr as an approach to suppressing APD altemans 

has been examined. iKr plays an important role in cardiac repolarization, increasing to 
a maximum during phase 3 of the action potential, as the channel recovers from 
inactivation, and then decreasing, as the electrical driving force decreases and as 
deactivation of the channel increases (Sanguinetti et al,, 'Two Components of Cardiac 
Delayed Rectifier K+ Current. Differential Sensitivity to Block by Class III 
Antiarrhythmic Agents," J G^^z Physiol 96:195-215 (1990); Rocchetti et al., "Rate 
Dependency of Delayed Rectifier Currents During the Guinea Pig Ventricular Action 
Potential,' VP/zy^za/ 534:721-732 (2001); Gintant, G., "Characterization and 
Functional Consequences of Delayed Rectifier Current Transient Repolarization," Am 
J Physiol 278:h806-H817 (2000); Gintant, G„ "Regional Differences in Ik Density in 
Canine Left Ventricle: Role of Iks in Electrical Heterogeneity," Am J Physiol 
268:H604-H613 (1995); and Clay et al., "A Quantitative Description of the E-4031- 
Sensitive Repolarization Current in Rabbit Ventricular Myocytes," Biophys J 
69:1830-1837 (1995), which are hereby incorporated by reference in their entirety). 
Because Ikf contributes minimally to the action potential plateau, increasing Ikf may 
have little or no effect on the Ca^"^ transient. Consequently, by increasing l^r it may 
be possible to suppress APD altemans without adversely affecting contractility. 
|0086] At present, however, there are no Ikf agonists available to test this 

hypothesis. To circumvent this problem, isolated ventricular myocytes were infected 
with an adenovirus expressing HERG, the gene that encodes the pore-forming domain 
of Ijcr, to increase HERG protein expression level and the corresponding Ikh current 
(Nuss et al., "Overexpression of a Human Potassium Channel Suppresses Cardiac 
Hyperexcitability in Rabbit Ventricular Myocytes," 7 C/m Invest 103:889-896 (1999); 
and Waldo et aL, "Effect of D-Sotalol on Mortality in Patients With Left Ventricular 
Dysfimction After Recent and Remote Myocardial Infarction, The SWORD 
Investigators. Survival with Oral D-Sotalol," Lancet 348:7-12 (1996), which are 
hereby incorporated by reference in their entirety). After verifying that Ikf had 
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increased post-infection, the myocytes were paced at rapid rates to determine whether 
increasing Iri by overexpressing HERG suppressed APD altemans. As shown in 
Figure 16, over expression of HERG markedly increased Ikf, as recorded during 
action potential clamp, and suppressed APD altemans during rapid pacing. Over 
expression of HERG did not, however, significantly reduce Ica? suggesting that this 
approach to suppression of altemans need not be accompanied by a reduction in 
contractility. 

[0087] The current paradigm regarding Ikf and VF is that blocking l^r is 

expected be anti-arrhythmic, secondary to prolongation of APD and refractoriness 
[Waldo et al., "Effect of D-Sotalol on Mortality in Patients With Left Ventricular 
Dysfunction After Recent and Remote Myocardial Infarction. The SWORD 
Investigators. Survival with Oral D-Sotalol. Lancet 348, 7-12 (1996); Hondeghem et 
aL, "Class III Antiarrhytlimic Agents Have a Lot of Potential But a Long Way to Go. 
Reduced Effectiveness and Dangeers of Reverse Use Dependence/' Circulation 
81 :686-690 (1 990), which are hereby incorporated by reference in their entirety]. 
However, this approach has not been successful in preventing VF and, moreover, is 
associated with pro-arrhythmia [Waldo et aL, "Effect of D-Sotalol on Mortality in 
Patients With Left Ventricular Dysfunction After Recent and Remote Myocardial 
Infarction, The SWORD Investigators. Survival with Oral D-Sotalol. Lancet 348, 7- 
12 (1996), which is hereby incorporated by referenced in its entirety]. Prolongation 
of APD, as reflected by an increased duration of the QT interval on the EGG, and an 
increase in the incidence of ventricular arrhythmias also are associated with the 
administration of numerous non-cardiac drugs that block iKr [Mohammad et aL, 
"Blockage of the HERG Human Cardiac K+ Channel by the Gastrointestinal 
Prokinetic Agent Cisapride," ylm. J. Physiol, 723:H2534-H2538 (1997); Roy et al, 
'TIERG, a Primary Human Ventricular Target of the Nonsedating Antihistamine 
Terfenadine," Circulation 94:815-823 (1996); Mitcheson et aL, "A Structural Basis 
For Drug-Induced Long QT Syndrome," Proc. Natl. Acad Sci. U.S.A. 97:12329- 
12333 (2000); Mitcheson et al., "A Structural Basis For Drug-Induced Long QT 
Syndrome," Proc. Natl. Acad Sci. U,S.A, 97:12329-12333 (2000), which are hereby 
incorporated by reference in their entirety]. Similarly, inherited loss-of-flinction 
mutations in Ikf are accompanied by a prolongation of the QT interval and by an 
increased risk of lethal ventricular tachyarrhythmias such as torsade de pointes 
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[Keating et aL, ''Molecular Genetic Insights Into Cardiovascular Disease," Science 
272:681-685 (1996); Curran et aL, ''A Molecular Basis For Cardiac Arrhythmia: 
HERG Mutations Cause Long QT Syndrome," Cell 80:795-803 (1995); Roden et al., 
"A Plethora of Mechanisms in the HERG-Related Long QT Syndrome. Genetics 
Meets Electrophysiology," Cardiovasc, Res, 44:242-246 (1999); Sanguinetti, et al, 
*'A Mechanistic Link Between An Inherited and An Acquired Cardiac Arrhythmia: 
HERG Encodes the kr Potassium Channel," Cell 81 :299-307 (1995), which are 
hereby incorporated by reference in their entirety]. 

[0088] The cardiac arrhythmias associated with inherited or drug-induced 

abnormalities of iRr are thought to be precipitated primarily by bradycardia-induced 
prolongation of repolarization [Keating et al., "Molecular Genetic Insights Into 
Cardiovascular Disease," Science 272:681-685 (1996); Curran et aL, "A Molecular 
Basis For Cardiac Arrh>thmia: HERG Mutations Cause Long QT Syndrome," Cell 
80:795-803 (1995), which are hereby incorporated by reference in their entirety]. The 
potential mechanisms by which Ikf rnay facilitate the induction of arrhythmias at slow 
heart rates have been studied extensively [Hondeghem et al., "Class III 
Antiarrhythmic Agents Have a Lot of Potential But a Long Way to Go. Reduced 
Effectiveness and Dangeers of Reverse Use Dependence/' Circulation 81:686-690 
(1990); Mitcheson et aL, "A Structural Basis For Drug-Induced Long QT Syndrome/' 
Proc, Natl Acad. Set U.S.A. 97:12329-12333 (2000); Roden et aL, "A Plethora of 
Mechanisms in the HERG-Related Long QT Syndrome. Genetics Meets 
Electrophysiology," Cardiovasc. Res, 44:242-246 (1999), which are hereby 
incorporated by reference in their entirety]. However, the contribution of Ikf to 
repolarization during tachycardia, which might also be important for the development 
of cardiac arrhythmias, has not been well characterized. The observation that 
reducing iKr increases the magnitude of APD alternans may provide an additional 
mechanism to account for the proarrhythmic effects of Ikf blockers, in that increased 
alternans magnitude would be expected to destabilize ventricular tachyarrhythmias, 
leading to the development of VF. Conversely, the observation that increasing l^r 
reduces the magnitude of APD alternans and the slope of the APD restitution relation 
provides a rationale for the development of a new class of compounds, iKr agonists, 
with the expectation that such compoimds may have anti-fibrillatory effects. 
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[00891 The latter hypothesis reUes, however, on the expectation that iKr can be 

increased sufficiently to reduce the slope of the restitution relation without shortening 
action potential duration to such an extent that contractility is impaired, secondary to a 
shortening of the plateau duration and attenuation of Ica- In addition, shortening of 
APD may lead to a reduction in the wavelength of reentry circuits (where wavelength 
= refractory period x conduction velocity). Reduction of the wavelength could, in 
turn, precipitate wavebreak by other mechanisms (e.g., so-called "head-tail" 
interactions, where a wavefront encounters a waveback and fragments) [Hund et al., 
"Dynamics of Action Potential Head-Tail Interaction During Reentry in Cardiac 
Tissue: Ionic Mechanisms," Am. J. Physiol. 279:H1869-H1879 (2000), which is 
hereby incorporated by reference in its entirety]. These and other potential drawbacks 
to the use of K-channel agonists to suppress VF remain to be evaluated critically. 
[0090] Until recently, therapy for the prevention of sudden cardiac death has 

been based on the presumption that frequent ventricular ectopy, in particular 
ventricular tachycardia, is a prelude to ventricular fibrillation [Myerburg et al., "Life- 
Threatening Ventricular Arrythmias: the Link Between Epidemiology and 
Pathophysilogy," In Cardiac Electrophysiology: From Cell to Bedside (AZipes, D.P. 
and Jalife, J. eds), pp, 723-731, WB Saunders Co., Philadelphia (1995), which is 
hereby incorporated by reference in its entirety]. Accordingly, drugs that suppress 
inducible or spontaneously occurring ventricular tachycardia are expected to prevent 
sudden death. However, recent large-scale clinical trials have indicated that classes of 
drugs that are effective for the suppression of ventricular tachycardia, e.g.. Class I and 
Class III antiarrh>1:hmic drugs, do not prevent sudden death and, in fact, may be 
proarrhythmic [Waldo et al., "Effect of D-Sotalol on Mortality in Patients With Left 
Ventricular Dysfunction After Recent and Remote Myocardial Infarction. The 
SWORD Investigators. Survival with Oral D-Sotalol. Lancet 348, 7-12 (1996); Echt 
et al., "Mrtality and Morbidity in Patients Receiving Encainide, Flecainide or Placebo: 
the Cardiac Arrhythmia Suppression Trial,'^ New Engl J, Med, 324:781-788 (1991), 
which are hereby incorporated by reference in their entirety]. In contrast, other 
classes of drugs that are not particularly effective for the suppression of most forms of 
ventricular tachycardia, such as beta-adrenergic receptor antagonists [Yusuf et al., 
"Overview of Results of Randomized Clinical Trials in Heart Disease: I. Treatments 
Following Myocardial Infarction," J, Am. Med. Assoc. 260:2088-2095 (1988), which 
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is hereby incorporated by reference in its entirety] and calcium channel antagonists 
[Held et al., 'Impact of Calcium Channel Blockers on Mortality," In Cardiovascular 
Pharmacology and Therapeutics (Singh, B.N. et al., eds), pp, 525-533, Churchill 
Livingston (1994), which is hereby incorporated by reference in its entirety], may 
reduce mortality from sudden death. 

I0091J If a causal relationship between the slope of the APD restitution 

relation and VF is confirmed, such a result could have significant implications for the 
pharmacological therapy of sudden cardiac death. Drugs that reduce the slope of the 
restitution relation would be expected to prevent the development of VF, but would 

not be expected to suppress ventricular tachycardia, if ventricular tachycardia is 
caused by some variant of spiral wave reentry. Conversely, drugs that do not reduce 
the slope of the restitution relation would not be expected to prevent VF, although 
they might suppress ventricular tachycardia, perhaps via a mechanism that does not 
involve alteration of restitution kinetics (e.g. slowing of conduction or prolongation of 
refractoriness). 

[0092] Given these expectations, the usefialness of existing agents for the 

prevention of VF could be reevaluated in light of their effects on the restitution 
relation and new drugs targeted against restitution could be developed. Naturally, 
such an effort would need to recognize that factors other than APD restitution 
contribute importantly to the induction and maintenance of VF. Furthermore, drugs 
intended to alter restitution selectively might have additional effects that could offset 
their intended effect. Nevertheless, judicious alteration of the APD restitution slope 
appears to be a promising approach to the treatment of VF and, as such, represents a 
potentially fruitful opportunity for drug development. 

[0093] Electrical restitution is a recent addition to factors that play a key role 

in the development of ventricular tachyarrhythmias. Experimental interventions that 
reduce the slope of the restitution relation have been effective in suppressing VF. 
These results encourage further investigation of altering restitution as a means of 
preventing sudden death, with the realization that clinically useful interventions that 
flatten restitution without having untoward effects have yet to be developed. 
[0094] Although the invention has been described in detail for the purposes of 

illustration, it is understood that such detail is solely for that purpose, and variations 
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can be made therein by those skilled in the art without departing from the spirit and 
scope of the invention as defined in the claims that follow. 



